AD-A118  085 
UNCLASSIFIED 


GEORGIA  INST  Of  TECH  ATLANTA  PRODUCTION  ANO  DISTRIBU— ETC  F/fi  12/1 
APPLICATIONS  OF  DUALITY  ANO  STOCHASTIC  DOMINANCE  IN  RELIABILITY— ETC  fU) 
JAN  82  M  PINE DO  N00014-80-K-0709 

PORC-82-1  NL 


Ft***  ■*  V 


Sf?2i 


/\ 


PDRC  Report  Series 
No.  82-1 
January,  1982 


APPLICATIONS  OF  DUALITY  AND  STOCHASTIC 
DOMINANCE  IN  RELIABILITY  THEORY 


Michael  Pinedo+ 


PDRC -82-1 


t  School  of  Industrial  and  Systems  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332 


wovvu. 

■fr.U- -  -  ~ 


This  research  was  supported  in  part  by  the  Office  of  Naval  Research  under 
Contract  No.  N00014-80-k-0709  and  in  part  by  the  National  Science  Foundation 
under  Grant  No.  ECS-8115344.  Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  U.  S.  Government. 


ABSTRACT 


J  In  this  paper  we  discuss  a  notion  which  is  the  dual  of  the  failure 
rate  and  is  defined  as  f(t)/F(t).  This  concept  has  come  up  several  times 
in  the  literature,  but  has  never  been  investigated  in  detail.  In  this 
paper  the  name  survival  rate  is  used.  We  show  that  in  many  models  of 
practical  interest  the  survival  rates  of  certain  random  variables  have  nice 
properties.  This  makes  it  possible  to  obtain  bounds  on  the  distribution 
functions  of  these  random  variables.  Moreover,  we  show  that  a  form  of 
stochastic  dominance  based  on  the  survival  rates  of  the  random  variables 
has  some  interesting  applications.  A 
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APPLICATIONS  OF  DUALITY  AND  STOCHASTIC  DOMINANCE 
IN  RELIABILITY  THEORY 

A  well-known  characteristic  of  the  distribution  of  the  lifetime  of  a 
component  is  its  failure  rate  r(t) ,  where  r(t)  la  defined  as  f (t)/(l-F(t>) . 
One  can  imagine  r(t)  to  be  the  probability  of  failure  of  a  component  during 
the  Interval  (t,  t+xj  given  that  the  component  is  still  alive  at  time  t, 
and  where  the  size  of  the  interval,  x,  is  very  small.  So 


r(t)  -  lim  — 
x  -*■  0  * 


1-tot)  -  F(t)  f(t) 

(1  -  F(t))  1  -  F(t) 


The  following  notion  can  be  considered  as  the  dual  of  the  failure  rate: 
The  probability  of  failure  of  a  component  during  the  interval  [t-x,  t) 
given  that  the  component  is  not  alive  at  time  t  and  where  again  the  size  of 
the  Interval,  x,  is  very  small.  As  this  is  the  conditional  probability  of 
the  component  surviving  to  time  t-x,  with  x  very  small,  we  will  use  the  term 
survival  rate  and  will. denote  it  by  s(t),  i.e. 


s(t) 


lim 
x  +  0 


1 

x 


F(t) 


This  concept  has  been  briefly  discussed  by  Barlow,  Marshall  and  Proshan 
(1963,  page  380)  and  by  Keilson  and  Sumita  (1980,  page  25).  However,  in 
neither  one  of  these  papers  it  was  given  a  name.  Barlow,  Marshall  and 
Proshan  (1963)  considered  distributions  with  f(t)/F(t)  decreasing  in  t, 
i.e.  distributions  with  a  Decreasing  Survival  Rate  (DSR) .  They  noted  that 
the  Decreasing  Survival  Bate  (DSR)  property,  just  like  the  Increasing  Failure 


Sate  (IFR)  property,  is  preserved  under  convolution.  This  was  shown  by 
observing  that  when  the  random  variable  X  is  IFR  the  random  variable  -X 
is  DSR. 

In  reliability  theory  distributions  with  monotone  failure  rates  are 
very  Important.  Increasing  Failure  Rate  distributions  have  some  well- 
known  and  useful  properties.  The  most  important  property  is  probably  the 
fact  that  IFR-ness  is  preserved  under  convolution.  A  second  property  is 
the  fact  that  the  lifetime  of  a  k  out  of  n  system  is  IFR  when  it  is  built 
up  of  components  with  independent  identical  IFR  distributions  (see  Barlow 
and  Proshan  (1975),  page  107).  It  is  well-known  that  F  being  IFR  is  equi¬ 
valent  to  log  (F(x))  being  concave  or  F  being  a  Polya  frequency  function 
of  order'2  (PF^) •  Barlow  et  al.  (1963)  also  observed  that  a  distribution 

F  is  DSR  if  and  only  if  log(F(x))  is  concave.  The  following  relationships 
hold  between  IFR  distributions,  DSR  distributions  and  distributions  of 
which  the- densities  are  PF^: 


(IFR) 

(DSR) 


In  order  to  get  some  intuitive  feeling  for  these  classes  of  distributions, 
observe  that  a  PF^  density  function  implies  that  the  density  function  is 
unlmodal,  while  IFR-ness  implies  that  there  are  no  downward  jumps  in  the 
density  function  and  DSR-ness  implies  that  there  are  no  upward  jumps  in 
the  density  function.  For  examples  of  PF^  densities,  see  problem  7  of 
page  79  of  Barlow  and  Proshan  (1975).  The  class  of  DSR  distributions  is 


fairly  large  as  it  comprises  the  class  of  DFR  distributions  as  well  as  the 


class  of  distributions  of  which  the  densities  are  PF^.  At  the  other  hand, 
one  can  show  easily  that  the  class  of  distributions  with  an  Inareasing 
Survival  Rate  is  empty.  Even  random  variables  on  a  finite  support  cannot 
have  an  Increasing  Survival  Rate. 

In  reliability  theory  a  distribution  F  is  called  Increasing  Failure 


Rate  on  the  Average  (IFRA)  if  (x) dx/ 1  is  increasing  in  t.  This  is  equi¬ 

valent  to  -(1/t)  log(l-F(t))  being  increasing  in  t.  A  dual  of  the  IFRA 


concept  can  only  be  defined  for  distributions  with  an  upper  bound  B.  This 
dual  property  then  implies  that  /s(x)dx/(0-t)  decreases  in  t.  It  is  clear 
that  this  concept  cannot  be  very  useful  in  practice  as  it  is  only  defined 


for  random  variables  on  a  finite  support. 

In  this  paper  we  will  give  some  examples  of  stochastic  models  where  the 
survival  rates  ?. '  ^oecif ic  random  variables  play  an  important  role  and  give 
valuable  additional  information  about  the  models.  In  Section  1  we  give 
examples  of  popular  models  where  the  random  variables  of  Interest  are  DSR. 

In  Section  2  we  discuss  a  form  of  stochastic  dominance  between  random  vari¬ 
ables  based  on  their  survival  rates  and  give  some  examples  where  this  type 
of  dominance  can  occur.  In  Section  3  we  present  bounds  on  the  distributions 
of  DSR  random  variables  based  on  the  knowledge  of  a  mean  and  a  percentile. 

In  the  last  Section  we  discuss  how  this  notion  can  be  used  in  the  analysis 
of  stochastic  models. 


1.  Examples  of  Random  Variables  with  a  Decreasing  Survival  Rate  (DSR). 

In  this  section  we  discuss  a  number  of  stochastic  models  in  which  the 
survival  rate  of  certain  random  variables  play  an  important  role.  The 


4 


proofs  are  In  some  cases  rather  easy  and  will  then  be  left  to  the  reader. 

(i)  The  Lifetime  of  a  Parallel  System.  Let  X^,,  i-l,.,.fn,  be  a  random 

variable  with  survival  rate  s  (t).  Let  Y«max(X.  , . . .  ,X  )  and  let  s(t)  denote 

i  in 

the  survival  rate  of  Y.  We  present  the  following  proposition  without  proof. 

Proposition  1.  The  survival  rate  of  the  lifetime  of  a  parallel  system  >*»n 

be  obtained  through  the  summation  of  the  survival  rates  of  the  individual 

n 

components,  i.e.  s(t)  -  T  s.(t)  for  all  t. 

1-1  1 

Recall  that  the  failure  rate  of  the  lifetime  of  a  series  system  can  be 
obtained  through  the  summation  of  the  failure  rates  of  the  individual  components. 

of  n  System.  Let 

if  component  1  is  functioning 
otherwise 

the  structure,  function  4  of  a  system  Be  defined 

If  the  system  is  functioning  under  x 
otherwise 

to  be  Independent  binary  random  variables  with 
we  define  the  reliability  function  h(p)  - 
*C$(X)-1)  -  E($X) ) .  For  any  coherent  structure  +  (not  necessarily  k  out  of 
n)  a  dual  structure  can  he  defined,  where  ■  1  -  4(l-x)  (see  Barlow 


(ii)  The  Lifetime  of  a  k  out 


and  let  5  »  (x^,...,xn).  Let 
as  follows 


If  X^,  i-l,...,n,  are  assumed 
P(X±-1)  -  p±  -  1-PCXj-O)  then 
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and  Proshan  (1975),  page  12).  Let  h^(p)  denote  the  reliability  function  of 
the  primal  system.  It  is  a  well-known  fact  that  F(t)  -  h^(F^(t),...,Fn(t)), 
where  F^(.t)  (»1-F^(t))  is  the  probability  of  component  i  reaching  age  t  and 
F(t)(*l-F(t))  is  the  probability  of  the  system  reaching  age  t.  It  can  also 


be  shown  easily  then  that  F(t)  -  h.  (F.  (t),...,F  (t)).  This  has  the  following 

♦d  1 

consequences:  If  for  the  primal  system  $  it  is  known  that  when  the  distribu¬ 
tions  F^,  1*1,..., n  are  IFR,  the  failure  rate  of  the  system  is  increasing  (and/ 
or  concave,  convex,  unlmodal),  then  the  distributions  F^  being  DSR  implies  that 
the  survival  rate  of  the  lifetime  of  the  dual  system  is  decreasing  (and/or 
concave,  convex,  unlmodal).  This  leads  us  to  the  following  proposition. 


Proposition  2.  The  lifetime  of  a  k  out  of  n  system,  built  up  of  compon¬ 
ents  with  DSR  lifetimes  that  are  i.l.d. ,  is  DSR.  . 


Proof;  This  follows  from  the  fact  that  IFR-ness  is  preserved  under  formation  ' 
of  a  k  out  of  n  system  with  i.l.d.  component  lifetimes  and  the  fact  that  the 
dual  of  a  k  out  of  n  system  is  an  n-k+1  out  of  n  system.  ■ 

(iii)  The  Number  of  Components  that  are  Uom  upon  System  Failure.  Consider 
a  system  where  the  lifetimes  of  the  components  are  i.l.d.  continuous  random 
variables.  Let  N  be  the  number  of  components  that  are  off  when  the  system 
goes  off.  El  Neweihi,  Proshan  and  Sethuraman  (1978)  and  Ross,  Shashahanl 
and  Weiss  (1980)  have  studied  the  properties  of  N.  The  random  variable  N 
is  discrete  and  the  failure  rate  is  defined  as 


k 

Ross  et  al.  (1980)  have  shown  that  N  is  1FRA,  i.e.  (1/k)  £  A.  is  nonde- 

J-l  3 

creasing  in  k.  Moreover  they  showed  that  if  the  system  has  nonoverlapping 


minimal  cut  sets,  i.e.  the  minimal  cut  sets  have  no  components  in  common, 
N  is  XFR,  which  implies  that  P(N-k|N  k  k)  is  nondecreasing  in  k. 


Proposition  3.  When  the  system  has  nonoverlapping  minimal  path  sets,  i.e. 
the  m-tn-tmal  path  sets  have  no  components  in  common,  N  is  DSR,  which  implies 
that  P(N«k|N  <  k)  is  noninareasing  in  k. 

Proof:  This  was  shown  by  Ross  et  al.  (1980)  through  a  duality  argument. 

The  reader  should  note  however  that  there  is  a  slight  error  in  this  paper; 
corollary  2  on  page  364  should  read  nonincreasing  Instead  of  nondecreasing.  B 

(iv)  The  Poisson  Shock  Models  A  system  is  subject  to  shocks  which  occur 
according  to  a  Poisson  process  with  rate  X.  The  i-th  shock  causes  a  random 
amount  of  damage  X  ,  where  the  X±,  1*1,2,...,  are  l.l.d.  with  distribution 
F.  The  system  falls  when  the  total  accumulated  damage  exceeds  a  given  thresh¬ 
old  x.  Let  Y  denote  the  lifetime  of  the  system  and  let  H  denote  the  distri¬ 
bution  of  Y.  This  model  has  been  thoroughly  analyzed  by  Esary,  Marshall 
and  Proshan  (1973).  They  proved  part  (a)  of  the  following  proposition. 

Proposition  4. 

(a)  When  the  distribution  F  is  DSR,  the  distribution  H  is  IFR. 

(b)  When  the  distribution  F  is  IFR,  the  distribution  H  is  DSR. 

Proof:  The  first  part  is  proven  in  Theorem  4.9  and  Theorem  3.2, of  Esary  et 
al.  (1973) .  The  proof  of  the  second  part  is  similar  with  regard  to  the  part 
corresponding  to  Theorem  4.9  and  easier  with  regard  to  the  part  corresponding 
to  Theorem  3.2.  ■ 

The  next  example  concerns  some  popular  queueing  models. 


(v)  The  Gl\G\l  and  Gl\D\c  queues .  Let  denote  the  service  time  of  cus¬ 
tomer  i,  T^  the  interarrival  time  between  the  i-th  and  the  (i+l)th  arrival 
and  D^  the  delay  of  customer  i  (the  time  customer  i  spends  in  queue).  It 
is  well-known  that  for  the  Gl|G|l  queue  under  the  FIFO  discipline  the  fol¬ 
lowing  recursive  relationship  holds: 

Dn+1  “  nax(0’  Dn  +  Sn  "  V* 

this  leads  us  to  the  following  proposition. 

Proposition  5.  If  1=1,2,...,  is  DSR  and  T^, ,  1=1,2,...,  is  IFR, 

D^,  1=1,2,...,  is  DSR. 

Proof:  This  follows  from  the  fact  that  -T^,  1=1,2,...,  is  DSR  (so 

D  +  S  -  T  is  DSR)  and  the  fact  that  the  maximum  of  two  independent 
n  n  n 

DSR  random  variables  is  again  DSR.  Clearly,  the  time  in  system  of  each 
customer  is  also  DSR.  ■ 

That  this  result  can  be  extended  to  the  GI|d|c  system  follows  from 
the  fact  that  under  the  FIFO  discipline  customer  c.l,  1=1,2,...  may  be 
assigned  to  station  c,  customer  c.l  -i-  1,  1=0,., 2,...  to  station  1,  customer 
c.l  +  2,  1=0, 1,2,...  to  station  2,  etc.  Each  station  may  then  be  viewed  as 
a  Gl|D|l  system. 

The  last  three  examples  concern  stochastic  scheduling  models. 

(vi)  The  Uakespan  of  Jobs  Subject  to  Series-Parallel  Precedence  Constraints. 
Consider  a  set  of  jobs  subject  to  precedence  constraints,  where  the  prece¬ 
dence  constraints  have  the  form  of  a  series-parallel  digraph.  For  a  dis¬ 
cussion  of  these  graphs  see  Lawler  (1978).  Lawler  defines  the  class  of 


series-parallel  graphs  recursively  as  follows: 

(a)  A  single  node  with  no  arcs,  e.g.  G  _  (.{!}*♦)*  series-parallel. 

(b)  If  -  (H^A^)  and  G2  «  0*2»A2^'  ''h®16  Ni  n  N2  “  $>  are  series- 

parallel,  then  G  -  G1  *  G2  “  u  u  A2  u  ^  x  m^)  is  also 

series-parallel;  G  is  said  to  he  a  series  composition  of  G^  and  G^. 

(c)  If  G2  -  and  G2  -  CN2>A2)  where  n  N2  "  *  are  sexfes“ 

parallel  than  G  ■  u  G2  »  CN^  u  ®2»A1  u  A2)  is  also  series-parallel; 

G  is  said  to  be  formed  by  the  parallel  composition  of  and  G2« 

(d)  Only  those  digraphs  which  can  be  obtained  by  a  finite  number  of  appli¬ 
cations  of  rules  a,  b  and  c  are  series-parallel. 

Proposition  6.  When  the  processing  times  of  the  jobs  are  Independent  ran¬ 
dom  variables  with  DSR  distributions,  and  the  capacity  to  process  jobs  Is 
unlimited,  the  distribution  of  the  make  span  of  the  project  is  DSR. 

Proof:  This  is  true  because  any  series  composition  of  sets  of  jobs  is  equi¬ 
valent  to  a  convolution  of  Independent  random  variables  and  any  parallel 
composition  results  In  a  random  processing  time  that  Is  the  maximum  of  two 
Independent  random  variables.  DSR-ness  is  preserved  under  both  operations.  ■ 

(vll)  The  Too  Machine  Plan  Shop  with  Zero  Intermediate  Storage.  In  this 
model  there  are  n  jobs  which  have  to  be  put  at  time  t"0  In  a  sequence  j^, 

j2 . j^.  Job  j^  is. processed  first  on  the  first  machine.  After  finishing 

its  processing  there  it  goes  to  the  second  machine  while  joh  J2  starts  its 
processing  on  the  first  machine,  etc.  When  job  j^  has  finished  ltd  proces¬ 
sing  on  die  first  machine  hefore  job  j^_^  has  finished  its  processing  on 
the  second  machine,  blocking  may  occur,  which  implies  that  job  j  ^  has  to 
remain  on  the  first  machine,  preventing  subsequent  jobs  of  being  processed 


on  the  first  machine. 


Proposition  7.  When  the  processing  times  of  all  the  jobs  on  the  two 
machines  are  DSR,  the  makespan,  i.e .  the  time  to  complete  all  the  jobs, 
is  DSR. 

Proof:  The  time  job  j^,  i“2 , . . . ,n,  occupies  the  first  machine  is  the  maxi¬ 
mum  of  the  processing  time  of  job  j^  on  machine  1  and  the  processing  time 
of  job  jj_^  on  machine  2.  The  makespan  is  the  convolution  of  the  times 
that  job  J^,  i"l, . . .  ,n,  occupy  machine  1  and  the  time  that  joh  J  occupies 
machine  2.  ■ 

Cviil)  The  21)0  Machine  Flaw  Shop  with  Infinite  Intermediate  Storage.  In 
this  model,  because  of  the  infinite  intermediate  storage  between  the  machines, 
blocking  does  not  occur.  When  job  j^  completes  its  processing  on  machine  1 
before  job  jj_^  completes  its  processing  on  machine  2,  job  j^  is  stored  in 
between  the  two  machines  and  machine  1  starts  processing  job  We  pre¬ 

sent  the  following  proposition  without  proof. 

Proposition  8.  When  the  processing  times  of  the  jobs  on  machine  1  are  deter¬ 
ministic,  not  necessarily  identical,  and  the  processing  times  on  machine  2 
are  DSR,  the  makespan  is  DSR. 

2.  Stochastic  Dominance  Based  on  Failure  Rates  and  Survival  Rates 

Lehman  (1955)  introduced  an  ordering  between  random  variables  which  he 
called  Monotone  Likelihood  Ratio  Ordering  (MLRO) .  Two  continuous  random 
variables  and  are  said  to  be  increasing  monotone  likelihood  ratio 


ordered  if  for  s  <  t 
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f2(t)  f2(s) 

Pinedo  and  Ross  (1980)  defined  two  random  variables  to  be  increasing 
failure  rate  ordered  if  for  all  s  <  t 

F  (t)  F  (s) 

— = - a  — - 

F2(t)  F2(s) 

They  showed  that  this  ordering  was  equivalent  to  r(t)  5  r2(t)  for  all  t, 
where  r^(t)  (r2(t))  is  the  failure  rate  of  X^(X2>.  Moreover  Pinedo  and 
Ross  (1980)  showed  that  Monotone  Likelihood  Ratio  Ordered  (MLRO)  ^Failure 
Rate  Ordered  (FRO)  "^Stochastic  Ordered  (SO). 

Keilson  and  Sumita  (1980)  Introduced  an  ordering  between  X^  and 
which  we  will  call  Survival  Rate  Ordering  (SRO) .  Two  continuous  random 
variables  X^  and  X2  are  said  to  be  increasing  survival  rate  ordered  if  for 
s  <  t 


Fl(t)  .  Fl(8) 


Using  the  Identity  F(x)  “  exp(  -J s(t)dt)  one  can  easily  show  that  this 
ordering  is  equivalent  to  s^(t)  £  s2(t)  for  all  t,  where  s^(t)  (s2(t>)  is 
the  survival  rate  of  X^(X2).  It  can  also  be  established  that  Monotone  Lik 
lihood  Ratio  Ordered  (MLRO)  ^Survival  Rate  Ordered  (SRO)  "^Stochastic 


Ordered  (SO) ,  which  results  in  the  following  chains  of  implications 
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MLRO 


^  SRO 


SO 


Two  random  variables  that  are  at  the  same  time  both  FRO  and  SRO  are  not 
necessarily  MLRO. 

In  the  remaining  part  of  this  section  Stochastic  Ordering  is  denoted 

by  £  ,  Survival  Rate  Ordering  by  £  ,  Failure  Rate  Ordering  by  £,  and 

so  sr  i 

Monotone  Likelihood  Ratio  Ordering  by  £  .  Two  examples  of  applications 

n 

of  these  orderings  are  now  given. 

(i)  Orderinga  Between  Lifetimes  of  k  Out  of  n  Systems.  Consider  a  k  out 
of  n  system  with  components  that  are  independent  and  Identically  distri¬ 
buted  with  distribution  F^.  Consider  a  second  k  out  of  n  system  with  com¬ 
ponents  that  are  independent  and  identically  distributed  with  distribution 

V 

Proposition  9.  If  F^  and  F^  are  Failure  (Survival)  Rate  Ordered,  Che  life¬ 
times  of  the  two  systems  are  Failure  (Survival)  Rate  Ordered  in  the  same 
sense  as  F^  and  F^. 

Proof:  We  will  prove  this  proposition  only  for  FRO  random  variables  as 
the  proof  for  SRO  random  variables  is  identical.  Let  X(t)  denote  the  fail¬ 
ure  rate  of  a  system  of  identical  components  with  failure  rates  r(t).  It 
is  well-known  that 


X(t)  -  r(t) 


p.h'(p)  | 

h(p)  |p  -  f(t) 
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where  h(p)  =  h(p. , . . . ,p)  is  the  reliability  function  of  the  system  as  defined 
in  example  (ii)  of  Section  1  (see  Barlow  and  Froshan  (1975),  page  109,  prob¬ 
lem  6).  It  can  be  shown  easily  that  p.h'(p)/h(p)  decreases  in  p.  The  pro¬ 
position  then  follows  for  the  FRO  case.  The  proof  of  the  proposition  for  the 
SRO  case  is  similar.  ■ 


(ii)  Optimal  System  Assembly.  Optimal  system  assembly  has  been  discussed 
several  times  in  the  literature,  see  Derman,  Lieberman  and  Ross  (1972,  1974). 
The  models  discussed  here  are  slightly  different.  Two  systems  are  considered: 
A  series  system  and  a  parallel  system,  both  with  n  component  locations.  In 


each  location  the  wear  on  the  component  installed  is  different.  Let  w^  denote 
the  wearout  rate  of  location  j ,  l.e.  w^  represents  the  expected  amount  of 
wear  sustained  by  a  component  in  location  j  per  unit  time.  There  are  n  non¬ 


identical  components  available.  The  random  variable  X^,  1=1,..., n,  with  dis¬ 
tribution  F^  represents  the  total  amount  of  wear  a  component  can  sustain  be¬ 
fore  failure.  Suppose  we  wish  to  minimize  the  probability  of  the  system  still 


being  alive  at  time  t.  Let  be  the  random  amount  of  wear  a  component  has  to 
undergo  in  location  j  during  a  period  t.  So  Y^  is  a  random  variable  with  mean 
Wj.t.  To  analyze  this  example  further  we  need  the  following  result  from 
Brown  and  Solomon  (1973). 


Proposition  20:  Assume  the  random  variables  and  X^  to  be  MLRO  with  E(X^)  > 
E(X£)  and  g(x,y)  is  a  real  valued  function  satisfying  g(y,x)  >  g(x,y)  for 
y  >  x.  Then 


«<W  *<W 


Proof:  See  Lemma  1  in  Brown  and  Solomon  (1973) 
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Consider  now  a  series  system  with,  n  component  locations.  There  are  n 
components  available  to  Install  in  these  n  locations  and  we  wish,  tp  maximi ze 
the  probability  of  the  system  reaching  age.  t. 


Proposition  11:  (a)  If  Y,  s  Y.  £ 

X  m  2  n 


5  Y  and  X.  X„  ,, ,  <  t 
m  n  X  X  2  i  |  2 


the  probability  of  the  series  system  reaching  age  t  Is  maximized  when  assign¬ 
ing  component  i  to  location  I,  for  i=l,...,n. 

00  “  Y1  Sar  Y2  \r  —  *sr  Yn  Bnd  ^  ^  X2  >m  —  ^  the  P^pbabiUty 

of  the  series  sytem  reaching  age  t  Is  maximized  when  assigning  component  I 
to  location  i,  for  1=1, . . . ,n. 


Proof:  The  proposition  can  be  shown  easily  through  an  adjacent  pairwise  inter¬ 
change  argument  and  proposition  IQ.  | 

Consider  a  parallel  system  with  n  component  locations.  Again  there  are 
n  components  available  to  Install  in  these,  n  locations  and  again  we  wish  tp 
maximize  the  probability  of  the  system  reaching  age  t.  • 

Proposition  12:  (a)  If  T.  s  Y,  *  ...  s  Y  and  X.  2  X,  2  ...X  the 

1  m  2  m  m  n  l  ar  2  ar  n 

probability  of  the  parallel  system  reaching  age  t  Is  maximized  when  assigning 
component  i  to  location  i,  for  1=1, ... ,n. 

If  Y1  *f  T2  *f  *'*  Yn  and  *i  -*2  Sm  *  *  *  \  the  probability  of 

the  parallel  system  reaching  age  t  is  maximized  when  assigning  component  i 
to  location  i,  for  1=1,..., n. 


Proof:  This  proposition  also  can  be  shown  easily  through  an  adjacent  pair 
vise  Interchange  argument  and  proposition  10.  ■ 


3.  Bounds  on  the  Distributions  of  DSR  Random  Variables 


In  reliability  theory  a  great  amount  of  research  has  been  done  on  the 
development  of  bounds  on  distributions  with  a  monotone  failure  rate.  Bar- 
low  and  Marshall  (1964,  1965,  1967)  developed  bounds  for  distributions  that 
are  IFR,  DFR,  IFRA,  distributions  with  PFj  densities  and  distributions  with 
densities  that  are  descreasing.  These  bounds  were  developed  based  on  a 
certain  amount  of  information  on  the  distributions.  The  following  possibi¬ 
lities  were  considered  with  regard  to  the  information  that  may  be  available; 

(a)  The  first  moment, 

(b)  the  first  and  second  moment, 

(c)  a  percentile  and 

(d)  conditions  on  the  failure  rate. 

In  this  section  we  explain  how  the  results  of  Barlow  and  Marshall 
can  be  extended  in  order  to  obtain  bounds  for  DSR  distributions.  Ve  will 
deal  first  with  nonnegative  random  variables  that  have  an  upper  bound,  l.e. 
random  variables  with  a  finite  support,  and  subsequently  with  random  vari¬ 
ables  that  have  no  upper  bound.  Observe  that  IFR  (DSR)  random  variables 
that  are  truncated  from  above  and  (or)  from  below  preserve  their  IFR-ness 
(DSR-ness) . 

Consider  a  nonnegative  random  variable  S  that  is  bounded  from  above 
by  U,  i.e.  P(S  <  0)  ■  1  and  has  a  given  mean  v.  Clearly,  when  S  is  DSR, 
the  random  variable  T  ■  U  -  S  is  IFR  with  mean  y  -  U  -  v.  Barlow  and  Marshall 
(1965)  developed  lower  bounds  for  the  survival  probabilities  of  the  random 
variable  T  under  the  assumption  that  y  is  known,  i.e. 
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P(T  >  t)  > 


e'th 

0 


for  t  <  u 
for  t  >  p. 


Transforming  these  bounds  into  bounds  on  the  survival  probabilities  of  the 
random  variable  S  gives  us  the  following  upper  bound: 


P(S  >  t)  < 


f 

j 1_e-(U-t)/(U-v) 


for  t  <  v 

for  t  >  v 


Barlow  and  Marshall  (1965)  developed  in  the  following  way  an  upper  bound  on 
the  survival  probabilities  of  T:  The  equation 


■  tj* w*dx 
0 


has  a  solution  v  If  and  only  If  t  >  p.  As  w  Is  a  function  of  t  we  write 
O  *-0 

w  (t).  This  solution  w  (t)  Is  unique  snd 
o  o 


!1  for  t  <  |i 

w. (t)  for  t  >  v 

o 

Barlow  and  Marshall  (1967)  developed  tables  for  the  values  of  w  (t).  The 

o 

following  lower  bound i  now  can  be  obtained  for  the  random  variable  S: 


P(S  >  t)  > 


1  -  wo(U-t) 
0 


for  t  <  v 


for  t  >  v 
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These  upper  and  lower  bounds  for  the  random  variable  S  are  depicted  in 
Figure  1. 


Figure  1 

One  of  the  results  above  can  be  extended  easily  for  random  variables 
that  have  no  upper  bound,  provided  besides  the  mean  also  a  percentile  £ 

P 

is  known,  i.e.  F(£p)  -  p,  where  p  is  close  to  1.  Then  it  can  be  shown  that 
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whether  or  not  random  variables  have  a  decreasing  survival  rate  as  this 
Implies  that  the  density  function  has  a  certain  "smoothness". 

Moreover  in  each  one  of  the  eight  examples  presented  in  Section  1 
it  appears  that  the  random  variables  that  are  DSR  also  may  be  IFRA.  For 
example  (i)  it  is  known  that  the  lifetime  of  a  parallel  system  is  IFRA,  pro¬ 
vided  the  lifetimes  of  the  individual  components  are  IFRA.  The  same  is  true 
for  any  coherent  system,  including  the  k  out  of  n  system  of  example  (ii) . 

Ross  et  al.  (1980)  have  shown  for  example  (ill)  that  the  number  of  components 
that  are  down  upon  system  failure  is  IFRA  for  an  arbitrary  system,  provided 
the  lifetimes  of  the  components  are  l.i.d.  For  the  Poission  Shock  model  of 
example  (iv)  it  is  known  that  for  an  arbitrary  F,  H  is  IFRA.  For  the  Gl|G|l 
queue  of  example  (v)  it  can  be  shown  easily  that  when  T^  is  DSR  and  S^  Is 
IFRA,  is  IFRA.  In  example  (vi),  when  the  processing  times  of  the  jobs 
have  IFRA  distributions,  the  distribution  of  the  makespan  of  the  project  is 
IFRA.  For  the  two  machine  flow  shop  without  intermediate  storage  of  example 
(vii)  the  same  is  true.  For  the  two  machine  flow  shop  with  infinite  inter¬ 
mediate  storage  of  example  (viii)  the  following  is  true:  When  the  proces¬ 
sing  times  of  the  jobs  on  machine  1  are  deterministic,  not  necessarily  iden¬ 
tical,  and  the  processing  times  on  machine  2  are  IFRA,  the  makespan  is  IFRA. 

Ve  would  like  to  finish  this  discussion  with  a  conjecture,  which  illus¬ 
trates  the  duality  between  IFR  and  DSR  random  variables  once  more:  Let 

X,,X0 . X  be  l.i.d.  IFR  (DSR)  random  variables,  then  Y  -  6(X,  ,X0, . . .  ,X  ) 

j.  i  n  x  i  n 

is  IFR  (DSR)  provided  the  function  $  is  Increasing  concave  (convex)  in  X^» 
M,...,n.  This  conjecture  can  be  easily  verified  for  n»l.  To  get  some 
additional  feeling  for  why  this  conjecture  may  be  true,  it  may  be  useful  to 
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j  recall  that  IFR-ness  (DSR-ness)  Is  preserved  under  convolution  and  that 

the  minimum  (maximum)  of  i.l.d.  1FR  (DSR)  random  variables  is  IFR  (DSR); 

|  the  minimum  (maximum)  is  an  increasing  concave  (convex)  function. 
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